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Abstract 



In the present study we calculate the energy values and the spatial distri- 
butions of the bound electronic states in some diffused quantum wells. The 
calculations are performed within the virtual crystal approximation, sp 3 s* 
spin dependent empirical tight-binding model and the surface Green function 
matching method. A good agreement is found between our results and ex- 
perimental data obtained for AlGaAs/GaAs quantum wells with thermally 
induced changes in the profile at the interfaces. Our calculations show that 
for diffusion lengths L D =20 -=- 100 A the transition (C3 - HHS) is not sen- 
sitive to the diffusion length, but the transitions (CI — HH1), (CI — LH1), 
(C2 — HH2) and (C2 — LH2) display large "blue shifts" as Lp increases. For 
diffusion lengths L D =0 4- 20 A the transitions (CI - HH1) and (CI - LHl) 
are less sensitive to the Ljj changes than the (C3 — HH3) transition. The 
observed dependence is explained in terms of the bound states spatial distri- 
butions. 

PACS numbers: 73.20.-r;7320.Dx 
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I. INTRODUCTION 



The Al x Ga,i_ x As / GaAs heterostructures have non-abrupt interfaces due to unwanted 
diffusion of Al and Ga across the heterojunctions. These compositionally graded interfaces 
change the electronic and optical properties of the quantum well structures [1-11]. A detailed 
study about the diffusion influence on the quantum well electronic states allows to control 
and use the new features these systems display. A comparison between calculated and 
measured transition energies permits to find diffusion coefficients for different temperatures 
[2,4,7,8]. Till now quantum wells with compositional grading at the interfaces were studied 
mainly using effective mass models and the envelope function approximation [1-4,6,8,9,11]. 
We believe that realistic tight-binding calculations will give a new information for these 
systems as was illustrated in [5] for an AlAs/GaAs quantum well with a linear concentration 
profile at the interfaces. 

In the present paper we study the influence of the diffusion length Lb on the optical 
transitions energies in an interdiffused Al^^Ga^As / GaAs quantum well with an as-grown 
well width of 100 A. This well was fabricated and experimentally studied in Ref.[l]. The 
energies of the transitions (CI — HH1), (CI — LH1) and (C2 — HH2) were measured at 
temperature 2 K for L D —0, 18, 20, 28, 48.5 and 66 A after annealing at 950 C for different 
time durations. 

II. MODEL AND METHOD 

The data of the system under study are taken from Ref . [1] . The width of the initial (as- 
grown) Al , 2 Gao,sAs / GaAs rectangular quantum well is 100 A. The Al composition across 
the well after the diffusion can be found from the expression: (See fig.l.) 

C(z) =C b + (C w - C 6 )/2 [erf (h - z)/L D + erf (h + z)/L D \ , (1) 

where C(z) is the Al concentration at a distance z from the well center; C w , the Al concen- 
tration in the well before the diffusion; Cf,, the Al concentration in the barrier before the 
diffusion; 2h, the width of the initial well; L D , the diffusion length. The Eq.(l) is widely 
used in the literature (see for instance Ref.[l]). 

We use the sp 3 s* spin dependent empirical tight binding model, the virtual crystal ap- 
proximation and the surface Green function matching technique [12-15]. For more details 
see [12] and references therein. The calculations were performed at the center of the two- 
dimensional Brillouin zone for the (100) growth direction. 

Each diffused structure was divided into three regions, namely, the external homogeneous 
barriers with negligible diffusion (the Al concentration in the planes of matching was 0.199) 
and the internal inhomogeneous (diffused) well region having a concentration profile as 
determined from Eq.(l). The size of the inhomogeneous slab was 167, 131, and 95 monolayers 
respectively for L£>=100, 70 and 35 A. The Green function of the external barriers was 
calculated from the transfer matrix in the usual way and the Green function of the diffused 
well region, by means of the algorithm established and used to study other heterogeneous 
structures [12,13,15]. 
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III. RESULTS AND DISCUSSION 



We calculated the energies of the first three bound electron states CI, C2, C3 (see fig.2a) 
and the first five hole states HH\, LH1, HH2, HH3 and LH2 (see fig. 2b) in the studied 
quantum well for different diffusion lengths Ld- The calculations in the interval Ld = 04 70 
A were performed through 5 A and in the interval L D = 70 4- fOO A - through fO A. We 
also found the bound states energies for the diffusion lengts Ld=0, 18, 20, 28, 48.5 and 66 
A. (The measurements in [1] were conducted with these diffusion lengths values.) As Ld 
increases the bound states depart further from the GaAs band edges. Qualitatively similar 
results have been reported in [5] for an Al As /GaAs quantum well with linearly graded 
interfaces. As a difference from [5], in our case the energies of the states C3, HH3 and LH2 
change several times less than the energies of the states CI, C2, HH1, HH2 and LH1. For 
instance, over the whole interval Ld = 4 70 A the state C3 rises only by 8 meV while over 
the same interval the state CI rises by 61 meV. The states HH3 and HH1 change by 19 
meV and 36 meV respectively over the same diffusion length interval. 

To understand the above mentioned behaviour we calculated the spatial distributions 
of all bound states for all diffusion lengths. The diffusion modifies the initial (as-grown) 
quantum well and leads to formation of two regions in the diffused well. In the as-grown well 
region (see fig. 1) the Al concentration increases and the energies of the bound electron (hole) 
states in the diffused quantum well increase (decrease). (The energy zero is at the Al As 
valence band top edge.) On the contrary, in the as-grown barrier region (see fig. 1) the Al 
concentration decreases and the energies of the bound states decrease (increase) for electrons 
(holes). All bound states are more concentrated in the as-grown well region and consequently 
their energies increase (decrease) for electrons (holes) as Ld increases. But, in general, a 
competition exists between two regions affected by the diffusion and a compensation of the 
energy shifts can occur. The magnitude of the energy shift for a given bound state depends 
on the relation between the probabilities of this state to occupy the diffused well part with 
increased Al concentration (the as-grown well region) and with decreased Al concentration 
(the as-grown barrier region). As an example, in order to illustrate the role of the above 
mentioned compensation, the spatial distributions of the electron states C3 and CI are 
shown in fig.3a and fig. 3b for the as-grown and a diffused (Ld=35A) quantum well. The 
energy of the state C3 changes only by 7 meV because this state spreads over both regions 
affected by the diffusion (within and outside the initial well, see fig. 3a) and the compensation 
occurs. The hole state HH3 has a similar spatial distribution and also shifts weakly. The 
energy of the state CI shifts with 21 meV due to its spatial localization in the central well 
region with increased Al concentration after the diffusion (see fig. 3b) and the compensation 
does not occur. The same argument is valid for the hole states HH1 and LH1. 

Theoretical curves of the transition energies dependence on the diffusion length Ld for 
the transitions (CI - HH1), (CI - LH1), (C2 - HH2), (C2 - LH2) and (C3 - HH3) are 
shown in fig.4a for Al . 2 Ga .$As / GaAs diffused quantum wells. There is a "blue shift" for 
all transitions with L D increasing as reported in [2,4] for similar systems. This shift is 27 
meV for (C3 - HH3) and 97 meV for (CI - HHl) as L D increases from to 70 A. The 
transitions (CI - LH1), (C2 - HH2) and (C2 - LH2) shift by 94 meV, 82 meV and 65 
meV respectively. 

The transition (C3 — HH3) is not sensitive to the diffusion length for Ld=20 4 100 A 
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because the energies of the states C3 and HH3 change weakly, see fig. 2. A relatively low 
sensitivity of the transition (C3 — HH3) to the diffusion length was measured in [2]. The 
transitions (CI - HHl), (CI - LHl) and (C2 - HH2) show considerable "blue shifts" in 
the interval L D =20 + 100 A because the energies of the states CI, C2, HHl, LHl and HH2 
change significantly, see fig. 2. 

For low diffusion length values (L D < 20 A) the transition (C3 — HH3) is sensitive to the 
Lo changes because the energies of the states C3 and HH3 depend strongly on the diffusion 
length, see fig. 2. The Al concentration changes only near the interfaces and the spatial 
distributions of the states C3 and HH3 do not allow a compensation. In this diffusion 
length interval the transitions (CI — HHl) and (CI — LHl) are less sensitive to the L D 
changes due to the energies of the states CI, HHl and LHl depend on L D weaker than the 
C3 and HH3 energies, see fig. 2. The states CI, HHl and LHl are mainly localized in the 
central well part where the Al concentration is not affected. 

In fig. 4b we compare our theoretical curves with the experimental curves from [1] for 
(CI- HHl), (CI -LHl) and (C2-HH2) transitions energies with L D =0, 18, 20, 28, 48.5 
and 66 A. The measurements in [1] have been conducted at temperature 2 K and temperature 
corrections of the theoretical values were not necessary. It is obvious that the agreement 
with the experimental data is quite satisfactory inspite of the fact that excitonic effects were 
not taken into account and fit was not done. All experimental curves lie about 10 -j- 45 
meV below the theoretical ones. The agreement between the theory and the experiment 
is better for the excited states than the ground states due to the fact that the exciton 
binding energies are higher for the ground states. For all transitions in fig. 4b two curves (a 
theoretical and the experimental) diverge with Lo increasing. For instance, for Ld=0 A the 
discrepances are 29 meV, 30 meV and 10 meV for the transitions (CI — HHl), (CI — LHl) 
and (C2 — HH2), respectively, and for Lp=66 A 45 meV, 41 meV and 29 meV respectively. 
A possible explanation of this behaviour could be an increase of the excitonic binding energy 
for higher Lp values due to the higher Al concentration within the initial well region. This 
means that the excitons are localized even in the case of Ld=QQ A. The spatial distributions 
of the bound states for these diffusion lengths support such an assumption. An exciton 
binding energy increase has been reported in [4] for double-barrier quantum wells as 
increases. 
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IV. CONCLUSION 



In the present paper we performed numerical calculations which showed that the empiri- 
cal tight-binding model combined with the surface Green function matching method and the 
algorithm for treating inhomogeneous finite slabs [12-15] works very well in case of quantum 
wells with graded interfaces. We found that the theoretical dependence of the transition 
energies on the diffusion lengths is in good agreement with the experiment. 

For given as-grown well width and diffusion length Lp, the sensitivity of the transitions to 
L D changes depends on the bound states spatial distributions over the diffused well regions. 
For small L D values (in our case below 20 A) the transition (C3 - HH3) is more sensitive 
to the L D changes than the transitions (CI — HH1) and (CI — LH1). For large L D values 
(in our case above 20 A) the transition (C3 — HH3) is not sensitive to the diffusion length, 
but the transitions (CI — HH1) and (CI — LH1) display large "blue shifts". 
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FIGURES 



FIG. 1. Al concentration profiles of quantum wells for diffusion lengths Ld=0, 25, 50 and 100 
A. The value of Ld=0 A corresponds to the initial as-grown rectangular quantum well of 100 A 
(35 monolayers). 

FIG. 2. (a) Electron energy levels (in eV) of the states CI, C2 and C3, (b) hole energy levels 
(in eV) of the states HHl, LH1, HH2, HH3 and LH2 for diffused quantum wells, as a function 
of the diffusion length Lo (in A). The value of Ld=0 A corresponds to the well width of 100 A 
(35 monolayers) of the initial as-grown rectangular quantum well. The energy zero is at the AlAs 
valence band top edge. 

FIG. 3. Spatial distributions of C3 (a) and CI (b) electron bound states over the as-grown 
(solid lines) and diffused (dotted lines) quantum wells. The diffusion length is Lo=3h A. The 
bound states energies are: C3 (as-grown well) - 2.240 eV, C3 (diffused well) - 2.247 eV; CI 
(as-grown well) - 2.118 eV, CI (diffused well) - 2.139 eV. 

FIG. 4. Transition energies (in eV) for diffused quantum wells, as a function of the diffusion 
length Ld (in A). The value of Le>=0 A corresponds to the well width of 100 A (35 monolayers) of 
the initial as-grown rectangular quantum well, (a) theory; transitions (CI — HHl), (CI — LH1), 
(C2 - HH2), (C2 - LH2) and (C3 - HH3). (b) theory (th.) and experiment (exp.); transitions 
(CI - HHl), (CI - LH1) and (C2 - HH2). 
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